Nucleon masses in strong magnetic fields 
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We compute the variation of the masses of the proton and the neutron induced by the presence of a strong 
external magnetic field. We discuss the choice of the wave function and different techniques how to apply the 
magnetic field. The results obtained from a 24 -16^ -12 and a 16^-8^ lattice are compared with a phenomenological 
approach. 
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1. INTRODUCTION 

We study the variation of the masses of the 
proton and the neutron induced by the presence 
of a strong external magnetic field. Our work 
is inspired by recent results of M. Bander and 
H.R. Rubinstein [1]. They consider the behavior 
of proton, neutron and electron in the presence of 
a constant magnetic field B. 

Their treatment assumes that for sufficiently 
small magnetic field the dependency of the nu- 
cleon masses on B is controlled by their anoma- 
lous magnetic moments which are known from ex- 
periment. The masses of proton and neutron can 
easily be computed using Landau's formula for a 
Dirac particle in a constant magnetic field [2]: 

E's=mp + ^il-9ps) (1) 



eB 

E" = mn + gns (2) 

2m„ 

s = ±1 denotes the spin, nip and m„ the nucleon 
masses for 5 = 0. The formulae hold for zero 
momentum, the lowest Landau level and eB > 0. 

Contributions to the masses which are even in 
the magnetic field such as \eB\ and |e5p are as- 
sumed absent in (2). Since such terms are not 
measurable by usual Zeeman experiments and 
cancel in the g-factor, the only way to estimate 
them is either experimentally (i.e. observing the 
conditions described in [1]) or by lattice simula- 
tions. Usually the presence of non-analytic terms 
of the type \eB\ is discarded in phenomenologi- 
cal formulae on physical grounds. However, such 



non-analytic terms do appear in the formulae for 
the point Dirac particle [2] and might survive in 
the low energy hadron masses too. 

Our simulations reported below are designed 
to put limits on these contributions. In fact, we 
confirm within the precision limits of our statis- 
tics the implicit hypothesis in [1] that such terms 
are absent. 

Ref [1] pointed out that if one inserts the exper- 
imentally known anomalous magnetic moments, 
gn = —1.9 and gp = 2.79, then for increas- 
ing B the lowest state of the neutron goes down 
steeper than the lowest state of the proton. Con- 
sequently, since m„ > nip there is a value of the 
magnetic field Bq where the masses of the two 
particles are equal. Bq is a very large field of 
the order of 10^'*T. Ref [1] also pointed out that 
for even larger fields (> IQ^^'T) a /3-decay of the 
proton may become possible. Fields of that order 
of magnitude may appear in certain astronomical 
phenomena. An unstable proton is certainly an 
interesting effect by itself. It also has implications 
for astrophysics (see [1] for references). 

However, the magnetic field required is very 
large, so that the assumption of pointlike parti- 
cles does not necessarily hold and the effect of the 
strong forces has to be taken into account. Since 
there is no analytic nonperturbative approach in 
QCD, it is worthwhile to treat the problem by 
means of lattice QCD. 

External magnetic fields have already been 
used as a tool to compute g-factors of proton and 
neutron [4]. In the present paper we study the 
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spin up/down states of proton and neutron and 
compare them which the phenomenological ap- 
proach of [1]. In particular we compute the slopes 
{energy I B) which is in a sense a harder task than 
computing the g-factors. For the g-factors it is 
sufficient to measure the mass difference of spin 
up and down states which is an observable that 
is more stable against systematic and statistical 
errors. In particular it is insensitive to terms even 
in the magnetic field, and to the QCD contribu- 
tions to the overall hadronic mass. 

2. THE METHOD 

We use a method similar to the one described 
in [4]. We work with a set of 30 quenched SU{?>) 
configurations separated by 500 sweeps with 2000 
sweeps thermalisation, [3 = 6.0. We use lattice 
sizes of 16^-8^ and 24 -16^ -12 and Wilson fermions 
at K = 0.1475 resp. K = 0.145. 

On the lattice a uniform magnetic field B can 
be introduced by multiplying a link Ufi(x) by a 
phase U^{x). For a field in z direction we set: 

B/~\ 



exp(—ieBa'^yLx) for x = — 1 
1 otherwise 
expi^ieBa"^ (x — xq)) 

where xq is an offset for the magnetic field to be 
discussed below. 

Consequently the plaquette in the x-y plane is 



Uy^{x) 



P^{x) 



exp(ie5a^(- 
for X 
expi^ieBa"^) 



otherwise 



Thus the magnetic fiux is homogeneous only if B 
is quantized: a'^eB = 2i:nl L^Ly . 

Unfortunately the lowest possible B is very 
large for reasonable sizes of the lattice. Due to 
the limited computer power available to us we de- 
cided to ignore the quantization. This introduces 
an inhomogenuity aX x = — l,y = Ly — 1. In 
order to limit its infiuence we put the source in 
the middle of the lattice and increased the lattice 
size in the x — t/-plane. 

We found that the masses depend on the choice 
of the offset xq in relation to the position of the 
source Xg . In order to get control about the pa- 
rameters we compared the 'fiat' case (all SU {"■!>) 



links set to one) with the predictions from the 
continuum. In particular we computed the sin- 
gle SU{Q) 'proton' and 'neutron' states and com- 
pared their dependence on B with the prediction 
given by the Landau formula [2] applied sepa- 
rately to each quark in the SU{Q) wave function. 
We expect that the proton spin up and the neu- 
tron spin down states should not depend on B 
while the states with opposite spin should go up 
for increasing B. 

If we compute these states for different choices 
for the offset xq we find a drastic difference be- 
tween Xq = Xg = LxI'i and x^) = 0. 




Figure 1. no QCD, W ■ 8^: 
«o = 0. = p|,n = p|,A = 



effective masses for 



Figs. 1 and 2 show the effective masses at t = 5 
(which is representative for all t). While the 
choice x^)=Xg reproduces the expected result the 
choice a;o = shows both proton states becoming 
heavier than the neutron states for sufficiently 
large B. This 'effect' can also be seen qualita- 
tively when QCD is switched on and may be con- 
founded with the physical effect predicted in [1]. 

In order to extract hadron masses from the cor- 
relation function the wave function of the particle 
of interest should have maximal overlap with the 
physical state. One way to achieve a better over- 
lap is certainly to use non-pointlike sources. Due 
to the interaction of the magnetic field configu- 
ration with the position of the source we have 
considered the simplest case of pointlike sources. 
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Figure 2. no QCD, 16^ • 8^: effective masses for 
for xq = Xg. O = □ = p|, A = rij, X = n|. 



We find that the overlap is much more sensitive 
to the choice of the wave function than it is for 
'non-magnetic' mass simulations. In fact we find 
the SU (6) wave function to have a better overlap 
than the (uj5Cu)d wave function [3]. 

3. RESULTS 

We computed the slopes (m(B) — m(0)) , 
the g-factors for proton and neutron and the ra- 
tio of their magnetic moments. The table shows 
the numerical results in comparison with the ex- 
perimental data. The data stem from a 16^ • 8^ 
lattice at K = 0.1475, a'^eB = 0.03, and from a 
24 • 16^ • 12 lattice at K = 0.145, a^eB = 0.02. 
They differ for the two lattice sizes, in particular 
for the spin up states. If we extrapolate to infinite 
volume by assuming a a = gq + l/Vai behavior 
we get results which are remarkably close to the 



experimental data. 

The ratios of the g-factors have been extracted 
from the expression [4] 
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(3) 



gn ' ' same but p ^ n 

Since this expression contains only ratios of 
masses we expect smaller systematical errors from 
the pointlike sources. Indeed ^(t) is constant 
within 2% for all timeslices t > 3. These results 
also hold for a^eB = 0.06 resp. 0.08. 

Our data for the slopes lead to the following 
conclusions: We find a dependency on the volume 
V . However, in the limit of infinite V the slopes 
give strong evidence that the linear phenomeno- 
logical approach of ref [1] holds in the strong field 
regime. 

In order to decide whether the proton becomes 
heavier than the neutron more precision is neces- 
sary. Further investigations thus should include: 

- a detailed study of the finite size effects 

- better statistics 

- smeared sources adapted to the problem 

- extrapolation to the chiral limit. 
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